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ABSTRACT 


Radar  polarimetry  is  a  relatively  recent  development  seeing  active 
research  only  in  the  last  few  decades.  The  phenomenon  that  optimal  (maximal 
power)  reflected  fields  may  exist  in  both  the  co-polarized  and  cross  polarized 
channels  of  the  receiving  radar  antenna  was  first  introduced  by  Kennaugh  and 
Huynen.  Current  research  efforts  focus  on  target  scattering  matrices  relating 
physical  attributes  of  the  target  that  bring  about  depolarization  of  the  scattered 
field  so  as  to  identify  the  various  scattering  mechanisms  and  to  enhance  target 
identification  and  detection. 

This  thesis  provides  a  comprehensive  survey  of  the  polarimetry  theories 
that  have  been  put  forth  by  various  researchers  to  characterize,  manipulate  and 
optimize  target  radar  returns  via  polarization  states.  These  polarimetry  theories 
include  matrix  representations  of  polarization,  the  various  vector  representations 
(Jones,  Sinclair,  Stokes,  Mueller  etc.)  and  their  usefulness.  This  thesis  also 
examines  the  Target  Decomposition  (TD)  theorem  that  seeks  to  decompose  the 
target  returns  into  individual  scattering  mechanisms  so  that  these  returns  can  be 
attributed  to  physical  characteristics  of  the  target.  Lastly,  the  topic  of  optimization 
of  polarization  states  for  maximizing  power  return  is  also  examined. 

Two  separate  models  are  implemented  in  Matlab  to  verify  and 
demonstrate  this  polarimetry  theory.  The  first  model  uses  TD  theorems  to 
simulate  foliage  clutter  and  study  its  effect  on  the  polarization  of  the  incident 
electric  field.  As  part  of  this  analysis,  a  (simulated)  static  dihedral  target  is 
introduced  and  its  effect  on  wave  polarization  is  also  simulated.  The  second 
model  was  constructed  to  study  optimization  of  polarization  states.  This  model 
takes  the  target  scattering  matrix  as  an  input  and  creates  all  possible  polarization 
states  of  incident  waves  to  interact  with  this  target.  It  then  graphs  the 
backscattered  power  in  the  co-polarized  and  cross-polarized  channels.  Both 
models  are  able  to  produce  the  expected  results  for  known  canonical  targets. 
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I.  INTRODUCTION 


Radar  polarimetry  is  a  promising  area  of  research  -  especially  when 
remote  sensing  applications  in  Synthetic  Aperture  Radar  (SAR)  are  considered. 
Traditionally,  radars  transmit  and  receive  waveforms  in  the  same  polarization 
state.  The  measured  data  include  the  amplitude,  frequency,  and  the  differential 
phase  of  the  received  wave  as  compared  to  the  transmitted  wave. 

However,  when  a  propagating  EM  wave  strikes  a  dielectric  or  a 
conducting  medium,  the  discontinuity  of  the  scatterer  will  generally  act  to 
depolarize  the  EM  wave,  possibly  resulting  in  a  completely  opposite  sense  of 
polarization  state  of  the  scattered  wave  to  which  the  radar  receiver  may  be  blind. 
All  information  (amplitude,  frequency  and  phase)  in  the  orthogonal  polarization 
state  of  the  back  scattered  wave  will  then  be  lost.  The  concept  of  depolarization 
can  be  illustrated  with  the  example  of  a  wire.  A  polarized  electric  field  that  is 
incident  on  the  wire  will  induce  currents  along  the  wire  and  these  currents  will  in 
turn  re-radiate  energy  that  is  polarized  in  the  direction  of  the  wire.  As  illustrated  in 
Figure  1-1,  the  scattered  wave  now  has  a  Ex  component  that  does  not  exist  in  the 
incident  wave.  If  the  incident  electric  field  is  orthogonal  to  the  direction  of  the 
wire,  one  would  not  expect  a  scattered  wave  since  no  currents  are  induced. 
Hence  the  polarization  of  the  scattered  electric  field  is  dependent  on  the 
geometry  or  orientation  of  the  wire. 


Figure  1-1  Example  of  a  wire  depolarizing  an  incident  electric  field  and 
scattered  wave  polarized  in  the  direction  of  wire,  with  components  in  the  x  and  y 

directions 
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A  polarimetric  radar  has  two  receiving  antennas  with  orthogonal 
polarization  states.  The  emitted  EM  wave  is  scattered  by  the  target  and 
measured  by  both  receivers  so  that  both  co-polarized  and  cross-polarized  signals 
are  determined.  These  measurements  provide  information  about  the  physical 
attributes  of  the  target.  Polarization  detection  may  also  allow  detection  of  some 
stealthy  targets  that  can  completely  depolarize  the  transmitted  wave,  resulting  in 
a  null  co-polarized  power  and  a  maximized  cross-polarized  power. 


A.  RECENT  RESEARCH  AND  DEVELOPMENTS 

The  basic  principles  of  radar  polarimetry  are  based  on  the  concept  of 
characteristic  polarization  states  first  introduced  by  Kennaugh  [1],  who 
demonstrated  that  there  exist  radar  polarization  states  for  which  the  radar 
receives  maximum  or  minimum  power  and  how  these  states  can  be  optimized. 
Hence  the  backscatter  power  of  the  target  depends  on  the  polarization  state  of 
the  incident  electric  fields.  Many  techniques  have  been  developed  to  determine 
the  optimal  polarization  state  and  involve  manipulation  of  the  target  scattering 
matrices.  The  results  are  usually  expressed  in  terms  of  the  power  received  by  the 
cross-polarized  and  the  co-polarized  channels  and  illustrate  the  importance  of  a 
full  polarimetric  data  set  for  deriving  target  information. 

Kennaugh’s  work  was  followed  by  Huynen’s  [2]  study  of  how  targets 
geometrical  properties  and  physical  structure  can  be  determined  by  their  back 
scattered  polarization  states.  Current  researchers,  including  Boerner  and  Cloude 
[1],  [2],  [4]  &  [5],  have  expanded  on  the  concept  of  Target  Decomposition 
theorems  that  serve  to  extract  information  about  the  target.  The  main  idea  is  to 
express  the  average  scattering  matrix  for  a  random  medium  as  a  sum  of 
independent  elements  and  to  associate  a  physical  scattering  mechanism  with 
each  of  these  independent  elements.  These  research  efforts  also  aim  to  integrate 
multi-polarization  data  into  multi-temporal,  multi-frequency  data  of  SAR 
interferometry  (POLINSAR)  so  as  to  achieve  greater  target  resolution,  clutter 
suppression  and  increased  detection  probability. 
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The  following  C-band  images  of  agricultural  fields  demonstrate  the 
dependence  of  the  radar’s  response  to  polarization  (Figure  1-2).  The  top  two 
images  are  like-polarized  (HH  on  left,  VV  on  right),  and  the  lower  left  image  is 
cross-polarized  (HV).  The  lower  right  image  is  the  result  of  displaying  these  three 
images  as  a  color  composite  (in  this  case,  HH  -  red,  VV  -  blue,  and  HV  -  green). 


Figure  1-2  Illustration  of  how  different  polarizations  (HH,  VV,  HV  &  color 
composite)  bring  out  different  features  in  an  agricultural  scene  [9].  (In  the  color 
plot,  HH  is  red,  VV  is  blue  and  HV  is  green) 

In  terms  of  actual  applications,  Table  1  gives  some  of  the  current 

developments  or  missions  employing  polarimetric  radars  in  remote  sensing. 
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Table  1-1  Polarimetrc  Radars 


Mission 

Application 

Polarization 

JPOLE 

Weather  radar  by  NSSL 

- 

SIR-C 

Remote  Sensing  on  Satellite 

HH,HV,VV,VH 

ALOS 

Remote  Sensing  on  Satellite 

HH,HV,VV,VH 

B.  THESIS  OUTLINE 

This  thesis  seeks  to  examine  and  verify  the  polarization  concepts  and 
models  put  forth  by  the  current  researchers.  The  author  will  make  use  of  the 
following  polarimetry  theories: 

1 .  Electromagnetic  wave  propagation  and  polarization; 

2.  Various  scattering  matrices  of  man-made  or  natural  objects  and 
their  relationship  to  the  physical  geometries  of  these  objects; 

3.  Target  decomposition  (TD)  theorems  and  their  basis  in  scattering 
physics.  TD  theorems  will  be  used  to  introduced  physical  variations 
into  scattering  matrices; 

4.  Optimization  of  polarization  states  for  various  target  scattering 
matrices  in  response  to  various  incident  electric  fields.  Polarization 
optimization  results  in  maximal/minimal  power  received  by  the 
receiving  antenna,  in  both  the  co-  and  cross-  polarized  channels 
and  this  is  unique  for  each  target. 

Theoretical  analyses  of  the  above  research  subjects  will  be  conducted  using 
Matlab  codes  to  investigate  the  behavior  of  various  targets.  The  analysis  will  be 
conducted  primarily  based  on  TD  theorems  and  polarization  optimization, 
respectively.  To  demonstrate  the  TD  theorems,  the  scenario  of  a  stationary 
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dihedral  scatterer  and  other  canonical  targets  hidden  in  foliage  clutter  will  be 
examined  and  polarimetry  will  be  employed  as  the  primary  tool  for  target 
detection.  To  examine  polarization  optimization,  an  analytical  technique  used 
currently  to  determine  optimal  polarization  states  will  be  implemented  in  Matlab. 
The  resulting  model  will  be  able  to  generate  the  optimal  polarization  states  for  the 
incident  electric  field  for  a  given  input  target  resulting  in  maximized/minimized 
power  in  the  co-polarized  and  cross-polarized  channels. 

Chapter  II  will  address  the  behavior  of  EM  wave  propagation  and 
interaction  with  surfaces.  Here  we  will  illustrate  the  idea  of  depolarization  from 
interaction  between  the  EM  wave  and  the  scatterer  due  to  the  physical  attributes 
of  the  scatterer. 

Chapter  III  illustrates  the  various  concepts  used  in  Polarimetry  Theory. 
These  include  tools  used  to  determine  the  polarization  states  of  scattered  EM 
waves  due  to  interaction  with  natural  or  man-made  scatterers  and  their 
theoretical  basis. 

Chapter  IV  describes  the  details  of  the  theoretical  analyses  conducted  and 
implemented  in  Matlab.  This  chapter  concludes  with  the  modeling  results  and 
their  analyses. 
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II.  ELECTRMAGNETISM  &  POLARIZATION 


A.  MAXWELL’S  EQUATIONS 

James  Clerk  Maxwell  (1831-1879)  put  together  the  laws  of 
electromagnetism  in  the  form  of  four  equations  based  on  theories  comprising 
previously  known  experimental  and  theoretical  results.  These  are  given  in  Table 
11-1.  Most  importantly,  he  introduced  the  concept  of  electromagnetic  waves  and 
the  physics  that  govern  their  propagation. 

Table  11-1  Maxwell’s  Equations  [16] 

Differential  Form  Integral  Form 


V.D  =  p,  ^ 

\>D-dS=  | 

)  L 

\pvdu 

) 

(2.1) 

VB  =  0  <j 

B-dS  = 

0 

(2.2) 

VxE  =  ~—  d 

dt  * 

\>E-dl  =  - 

—  \BdS 
dt{ 

(2.3) 

VxH  =  J  +  —  d 

dt  ^ 

:  dD 

( J  +  ~)  ‘  dS 
,  dt 

(2.4) 

Assuming  a  time  factor  of  e  jtM  and  noting  that5  =  pH ,  J  =  <jE  ;  and  taking  the 
curl  of  both  sides  of  the  differential  form  of  (2.3)  gives 

VxVx£  =  -ja>juS7xH  (2.5) 

V(V  •  is)  -  V2E  =  -jcoju(<r  +  jcoe)E  (2.6) 

V2E-y2E  =  0  (2.7) 

where 


y2  =  j  copier  +  jcos) 


(2.8) 


7 


and  y  is  called  the  propagation  constant.  Note  that  the  e  Ja*  disappears  because 
it  is  associated  with  every  term  and  therefore  factors  out,  resulting  in  time 
independent  equations.  By  a  similar  process,  it  can  be  shown  for  the  magnetic 
field,  H,  that 

V2H-y2H  =  0  (2.9) 

Equations  (2.7)  and  (2.9)  are  known  as  the  Helmholtz  equations.  The  solution  to 
the  first  equation  gives  the  electric  field  propagation  equation.  We  can  assume 
that  the  wave  propagates  along  the  +z  direction  and  that  E  has  only  an  x- 
component,  then 

E  =  Ex(z)  (2.10) 

Substituting  into  equation  (2.7)  gives 


(V2-y2)Ex(z)  =  0 

T2-r2]z,(z)  =  o 

az 


(2.11) 


Finally,  realizing  that  this  is  a  linear  homogeneous  differential  equation,  we  can 
insert  the  time  factor  and  obtain  the  final  solution  as  follows: 

E(z,t)  =  E0e~az  cos(cot-kz)  (2.12) 

where  E0  is  a  scalar  constant,  a  is  the  attenuation  constant,  co  is  the  angular 

frequency  and  k  is  the  phase  constant  or  the  wave  constant.  Equation  (2.12)  is 
the  electric  field  vector  that  defines  the  propagation  of  the  electric  field. 

The  magnetic  field  H  is  related  to  the  E-field  and  follows  by  a  similar 
argument.  In  a  different  (but  equivalent  form),  we  write 

H(z,t)  =  Re(H0e~azej(O),-kz) )  (2.13) 


H  =  ^-e  “zcos (cot-kz  +  S) 

\v\ 


(2.14) 
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where  rj  is  the  intrinsic  impedance  and  may  be  complex.  This  complex  nature  of 
?i  introduces  a  phase  difference  8  between  E  and  H  fields  (note  that  this  is  a 
phase  difference,  assuming  E-field  has  zero  absolute  phase).  In  free  space,  the 
E-field,  H  field  and  the  direction  of  wave  propagation  are  all  orthogonal  to  each 
other  and  can  be  represented  as  in  Figure  11-1 . 

E 


Direction  of  wave 
propagation,  z 

Figure  11-1  Electromagnetic  Wave  Propagation  in  Free  Space 

B.  WAVE  PROPAGATION  AND  POLARIZATION 

Having  defined  the  electric  and  magnetic  field  vector,  we  now  turn  our 
attention  to  the  interaction  of  EM-waves  with  various  media.  The  relative 
amplitude  of  the  incident,  reflected  and  transmitted  waves  depends  on  the 
constitutive  parameters  of  the  two  media  involved,  i.e.  the  permeability^, 
permittivity,  //  and  the  conductivity,  a .  Objects  that  tend  towards  conductors 
reflect  most  of  the  incident  wave  energy,  while  objects  that  tend  towards  lossy 
dielectrics  absorb  larger  portions  of  the  incident  energy. 

In  addition  to  energy  reflection  and  transmission,  the  waves  can  become 
depolarized  so  that  the  polarization  state  of  the  incident  wave  differs  from  that  of 
the  scattered  wave.  Hence  the  resultant  state  of  the  scattered  wave  depends  on 
the  polarization  of  the  incident  wave  as  well  as  the  parameters  of  the  reflecting 
medium.  We  discuss  the  elements  of  wave  propagation  in  the  following  sections. 
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1.  Polarization 

Consider  an  EM  wave  propagating  out  of  the  page  and  situated  at  the 
origin  of  an  x-y  axis  system,  as  shown  in  Figure  11-2.  The  E-field  is  represented 
by  the  vector  shown  and  varies  continuously  in  magnitude  and  changes  direction 
every  half-period  according  to  the  sinusoidal  behavior  described  in  Section  II.A. 
The  components  of  the  E-field  in  the  x-y  directions  are: 


^5 

II 

+ 

(2.15) 

Ex  (z,  t)  =  E0x  cos  {cot  - kz  +  8x ) 

(2.16) 

Ey  (z,t)  =  E0y  cos  (cot  -kz  +  5y) 

(2.17) 

y 


X 


Figure  11-2  Representation  of  instantaneous  E-vector 
Examining  the  two  components  Ex  and  Ey,  shows  that  the  sum  of  two  sinusoidal 
waves  traces  an  ellipse.  Their  mathematical  relationship  is  described  in  the 
following  equation  of  an  ellipse: 

c  F  ^  EE  ( E  ^ 

-2  *  *  cos 8  +  —  =  sin2  8  (2.18) 

J  E0xE0 y  J 

where  5  =  5y  -  Sx ,  is  the  phase  difference  between  Ex  and  Ey. 

Depending  on  the  parameters  of  the  ellipse  equation,  the  E-vector  may  be 
reduced  to  the  special  cases  of  linearly  polarized  (as  in  Fig  11-2),  circularly 
polarized  and  the  most  general  case  being  elliptically  polarized.  Linear 
polarization  occurs  when  the  two  orthogonal  components  Ex  and  Ey  are  in  phase 
and  so  the  direction  of  the  linear  polarization  depends  on  the  relative  amplitude 
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of  the  two  components.  Circular  polarization  occurs  when  the  two  orthogonal 
components  have  the  same  amplitude  but  are  exactly  ninety  degrees  out  of 
phase.  All  other  cases  (with  the  two  components  not  in  phase  and  not  having  the 
same  amplitude)  are  elliptical  polarizations. 


Figure  11-3  Illustration  of  (i)  Linear,  (ii)  Circular  and  (iii)  Elliptical  Polarization[19] 
The  polarization  ellipse  is  a  common  representation  of  the  state  of 
polarization.  The  tilt  angle,  cp  of  the  ellipse  is  obtained  from 


tan  2  <p  = 


2£oAyCQS^ 

E2  -E2 

ox  oy 


(2.19) 


The  ellipticity  angle,  z  is  simply  derived  from  the  major  and  minor  axes  of  the 
ellipse  as  follows, 


Minor  axis 

tan  t  = - 

Major  axis 


(2.20) 
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From  Mott  [11],  the  wave  can  also  be  characterized  by  more  tractable 
mathematical  terms  using  the  polarization  ratio,  which  is  a  complex  number, 
shown  here  with  its  relation  to  the  polarization  ellipse  parameters  [1] 


P 


A )y  is  .  js  tan  (p  +  tan  r 

— —  •  e  =  tan  <peJ  = - 1 - 

E0x  1  -  j  tan  (p  tan  r 


(2.21) 


For  a  circularly  polarized  wave,  the  EM  wave  is  a  sum  of  equal  amplitude 
left-hand  polarized  and  right-hand  polarized  components,  EL  and  ER.  The 
polarization,  after  some  manipulation,  is  then  [11] 


q 


1  -jP 
1  +  jP 


(2.22) 


Using  a  modified  polarization  ratio,  p  =  jP ,  we  then  obtain, 


q  = 


i-p 
i  +  p 


p  = 


i -q 

1  +  C[ 


(2.23) 


The  various  representations  provide  varying  degrees  of  ease  in 

representing  the  polarization  of  a  propagating  wave  and  are  also  utilized  for  more 
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complex  types  of  representations,  e.g.  the  Poincare  sphere.  A  more  in-depth 
treatment  of  polarimetry  theory  and  its  applications  to  radar  detection  will  be 
given  in  chapter  III. 

2.  Wave  Reflections  from  Smooth  Flat  Surfaces 

Using  a  plane  incident  wave  to  illustrate  smooth  surface  reflection,  the 
reflected  wave  is  also  a  plane  wave  and  this  is  known  as  specular  reflection.  The 
electric  field  of  the  reflected  wave  will  depend  on  the  polarization  of  the  incident 
wave.  It  is  convenient  to  assume  that  the  incident  wave  has  two  components,  as 
before.  One  component,  EiV,  is  the  “vertically”  polarized  component  and  the 
other,  EiH  is  the  “horizontally”  polarized  component.  This  is  illustrated  as  follows. 


Figure  11-5  Plane  wave  Reflection  from  a  Smooth  surface  (parallel  polarization) 


The  reflected  components  are  related  as  follows  [13] 

ErV  =  EVEiV 


ErH=THEiH 


(2.24) 


where  r  is  the  reflection  coefficient  and  depends  on  polarization  of  the  incident 
wave.  The  polarization  of  the  incident  wave  will  have  components  parallel  or 
perpendicular  to  the  plane  of  incidence  and  it  is  convenient  to  express  the 
associated  reflection  coefficients  separately  as  r;/  and  T±.  These  coefficients 
are  defined  as 
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(2.25) 


_Hr  _ri2  cos  6t  -  r/x  cos  0i 

n  —  — 

Ht  rj2  cos  6t  +  t/x  cos  6t 

r  _Er  _  t]2  cos  6t  -  //,  cos  0t 
1  Ej  r]2  cos  0j  +  r/x  cos  Ot 

where  7  is  the  intrinsic  impedance.  This  expression  shows  clearly  the 
dependence  of  wave  reflection  on  material  properties  and  wave  polarization. 


3.  Composite  Reflections 

For  smooth  but  irregular  surfaces,  each  part  of  the  incident  wave 
corresponds  to  reflections  in  different  directions.  As  a  result  a  divergence  factor 
D  j  suggested  by  Kerr  (1964)  [13],  needs  to  be  included  in  the  derivation  of  the 
reflected  E-field. 


In  the  case  of  a  rough  surface,  the  irregularities  can  be  thought  of  as 
fluctuations  about  an  “average”  flat  value  and  such  surfaces  can  still  reflect 
waves  in  a  specular  manner.  Rough  surfaces  can  be  modeled  by  assuming  a 
Gaussian  variation  of  surface  height  and  have  a  roughness  (loss  coefficient)  of 
[13], 


Ps  =  eXP 


( nh  ^ 

2 

-8 

rms 

sin2  6 

K  A  j 

(2.26) 


where  hrms  is  the  standard  deviation  of  the  height  variations  about  the  average 
flat  surface  and  6  is  the  incident  angle. 

By  combining  results  of  the  preceding  sections,  we  may  determine 
magnitude  and  phases  of  the  reflected  fields  when  reflections  occur  from  a  rough 
spherical  surface,  as  [13]: 


ErV=TvDpsEiv 

ErH=THDpsEiH 


(2.27) 
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C.  RADAR  CROSS  SECTION  (RCS)  a 

The  radar  cross  section  is  the  portion  of  the  back-scattered  power 
corresponding  to  the  particular  polarization  of  the  receiving  antenna.  Thus, 


a  =  lim  4  7tR2 

R — ^oo 


power  per  unit  area  in  scattered 

wave  at  receiving  antenna 
<  > 
which  is  in  the  polarization  of 

receiving  antenna 

f  power  per  unit  area  in 

\wave  incident  on  target 


This  definition  implies  that  only  the  portion  of  the  scattered  wave  at  the  receiving 
antenna  that  has  the  same  polarization  of  the  antenna  is  used  in  defining  a . 
Since  the  receiving  antenna  may  have  any  specified  polarization,  the  definition  of 
cr  of  any  arbitrarily  polarized  scattered  wave  should  be  divided  into  two  parts  - 
one  to  which  the  radar  responds  and  the  other  that  is  orthogonal  to  the  first. 

We  define  the  E-field  scattered  by  the  target  and  returning  to  the  receiving 
antenna  (in  cylindrical  coordinates)  by  the  sum  of  components 


E S  —  EgO  +  E Jj)  — 


(2.28) 


Figure  11-6  gives  the  geometrical  relation  between  the  target  and  the  receiver  and 
the  definitions  of  the  various  parameters. 


Figure  11-6  Scattering  Geometry 
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Using  matrix  notation  and  defining  gas  the  polarization  ratio  [13],  which 
carries  the  polarization  information  of  the  EM  wave  and  is  analogous  to  P  for  the 
case  of  plane  waves, 


the  E-field  becomes, 


(2.29) 


(2.30) 


Next  we  define  the  scattered  E-field,  Es  as  a  composition  of  polarized 
states,^,  having  an  arbitrary  but  specified  polarization  and,  EB  that  is 
orthogonal  to  EA .  These  components  are  represented  as  [13] 


Ea  -  A@  +  AxQj>  -  A x 


1 

0, 


EB  —  A^O  +  Az(j) 


(2.31) 


E  -  F  +  E 
n,s  —  r,  A  ^  n  B 


where  Ax,  A  ar|d  4  are  parameters  to  be  defined  so  that  the  above  is  true,  and 
4*0  is  assumed.  QA  is  the  polarization  ratio  based  on  the  polarization  of  the 
receiving  antenna.  Since  EA  is  orthogonal  to  EB ,  the  following  condition  must  be 
satisfied: 


Ea'Eb  =  Ax[  1  qa\ 


=ax(a;+qaa;)=o 


(2.32) 


Using  equation  (2.32)  and  (2.31)  to  solve  for 4,  4  ar|d  4>  we  obtain  the 
following  expressions, 
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i  M±gg) 

>+laf 


,  -^gxg-g,) 
2  —  2 
l  +  l&l 


(2.33) 


,  UQzM 

3  H&f 

This  result  shows  that  any  arbitrarily  backscattered  wave  can  be 
decomposed  into  the  sum  of  2  waves,  EA  and  its  orthogonal  component  EB ,  to 
which  the  receiving  antenna  can  respond.  The  power  density  in  the  wave 
component  EA,  which  will  be  received  by  the  antenna  is  then, 


Power  density  in  EA  = 


\ea\2  e8\  i+q:q 


2/7  2/7(1  + |a  |  ) 


(2.34) 


We  can  express  a  in  terms  of  the  polarization  ratios  and  as  a  function  of 
the  total  scattered  RCS,  cr  . 


Es 

|2  I 

L  =  lim  4/rR2  - 

E  7 

^ 9 

E, 

R2  — >oo 

E,2 

(i+|gf> 


g  =  lim  4 7rR2 

s  i?2->  00 


|1  +  ae|' 


R2  — >oo 


\E\ 


Rj  — >oo 


\E,  I  (H& 


(2.35) 


Finally,  we  obtain, 


cr  =  cr. 


Ppo/  — 


Ppol 


^  +  QaQ\ 


a+|a|  )a+|e| ) 


(2.36) 
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ppol  is  known  as  the  polarization  efficiency  or  polarization  match  factor  [11], 

having  a  range  of  0<  ppol  <1 .  The  polarization  match  factor  shows  how  well  a 

receiving  antenna  is  matched  in  polarization  to  an  incoming  wave.  It  can  be  used 
to  study  scenarios  whereby  the  polarizations  of  the  transmitting  and  the  receiving 
antennas  are  different. 

D.  THE  RADAR  SYSTEM 

This  section  will  describe  the  operation  of  a  radar  system  using  some  of 
the  aforementioned  theories. 

1 .  The  Radar  Signal  Waveform 

Similar  to  the  EM  wave,  the  radar  signal  waveform  takes  the  form  of  a 
time-varying  sinusoidal  wave  that  may  be  represented  as  [10], 

s{t)  =  a(r)cos[<sy +  #(r)  +  $y  (2.37) 

where  a  gives  the  amplitude,  6{t)  is  the  phase  term,  is  a  arbitrary  phase 
angle  and  a>0  is  the  carrier  frequency.  Usually,  the  carrier  frequency  is  removed 

before  signal  processing.  Hence,  amplitude,  phase  and  frequency  modulation  are 
directly  related  to  these  parameters.  The  (modulus  of  the)  Fourier  transform  of 
the  radar  signal  gives  the  spectrum  of  the  signal. 

The  signal  waveform  is  often  defined  over  the  time  interval  for  which  the 
radar  signal  is  pulsed.  A  pulse  train  of  waveforms  is  illustrated  in  Figure  11-7,  with 
the  pulse  width  and  the  pulse  repetition  interval  (PRI)  defined  graphically. 


Figure  11-7  Pulse  train 
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The  average  peak  transmitted  power  is  given  by, 


The  average  power  is  then, 


(2.38) 


where  t  is  the  pulse  duration  and  T  is  the  integration  time. 


(2.39) 


2.  Radar  Equation  [10] 

The  power  transmitted  by  the  radar  antenna  is  first  assumed  to  be 
radiated  spherically,  and  the  power  density  at  a  distance  R  from  the  radar  is  then, 


Pt 

4kR 2 


(2.40) 


where  Pt  is  the  power  output  by  the  transmitting  antenna. 


Taking  into  consideration  a  directional  antenna  in  order  to  increase  power 
density  in  a  certain  direction,  the  power  density  is  multiplied  by  a  directive  gain,  G 
that  is  usually  a  scalar  function  of  the  direction  angles. 


Pfi_ 

4kR2 


(2.41) 


The  radiated  power  is  then  incident  on  the  target  and  the  reflected  amount 
is  determined  by  the  RCS,  <j  .  The  reflected  power  is  re-radiated  spherically  from 
the  target  and  a  behaves  like  an  “antenna  gain”. 


Pfi  o' 

4 kR2  4k R2 


(2.42) 


The  reflected  power  is  then  intercepted  by  the  radar  antenna,  with  effective  area 


A.  = 


GA2 

4k 


and  so  the  final  power  received  is, 


PtGo-  _  PtG2 A2 a 

(4k R2)2  e  ~  (4kY  R4 


(2.43) 
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In  practical  situations,  the  returned  signals  received  by  the  radar  will  be 
corrupted  by  thermal  noise, 

Nt=kTB  (2.44) 

where  k  is  the  Boltzman’s  constant,  T  is  the  effective  temperature  and  B  is  the 
bandwidth  (equal  to  1  It).  In  order  to  derive  the  Signal-to-Noise-Ratio  (SNR),  we 
let  the  power  received  be  the  minimum  detectable  signal  by  the  radar,  Smin .  The 

quantities  Smm  ,  Ni  and  SNR  are  related  by  a  radar  figure  of  merit,  F ,  the  noise 
figure. 


S . =  kTBF  (SNR) 


Finally,  the  radar  equation  describing  the  SNR  becomes, 


(SNR)  = 


(■ AnfkTBFLR 4 


(2.45) 

(2.46) 


where  L  denotes  the  radar  losses  [10]. 
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III.  POLARIMETRY  THEORY 


A.  JONES  VECTOR 

The  Jones  vector  is  introduced  as  a  matrix  representation  of  the 
polarization  state  of  the  electric  field.  The  Jones  vector  that  represents  a 
polarization  state  can  be  decomposed  into  any  orthogonal  polarization  basis. 

Using  the  definitions  in  equation  (2.15),  a  propagating  EM  wave  E(z)  is 
represented  as  a  sum  of  its  components  in  the  x  and  y  directions.  Because  these 
components  oscillate  with  the  same  frequency  for  the  monochromatic  case,  we 
can  suppress  the  temporal  term,  exp (jcot )  and  the  equation  becomes, 


E(z)  =  e-jkz 


Ex 

E„ 


(3.1) 


Furthermore,  as  the  wave  is  planar,  the  magnitude  of  the  electric  field  E(z) 
is  the  same  at  every  point  of  its  propagation.  Hence  we  may  set  z=0  for 
convenience  and  the  equation  is  further  reduced  to, 


m=£0 


(3.2) 


The  vector  E(0)  is  called  the  Jones  vector  of  the  wave.  The  amplitude  and 
phase  of  the  complex  components  of  the  electric  field  are  completely  defined  by 
this  vector. 


1.  Ellipticity  Definitions 

Once  again,  we  can  define  the  phase  difference  between  the  components 
of  the  field  as  8  =  8  -  8 .  Then  from  the  Jones  vector,  we  can  determine  the 

y  x 

information  necessary  to  construct  the  polarization  ellipse  [8]. 
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The  tilt  <P  is 


2  I  E  I 

tan  2cp  = - - — £ - -  cos  8 

T  I  Z7  |2  I  77  |2 


(3-3) 


The  ellipticityr  is 


sin  It  = 


2  £  £ 


71  sinJ 


2  2 

^1  +  E, 


(3-4) 


The  magnitude  of  the  wave  is 

e „ = Vi  T  i2  + 1 T 


(3.5) 


Conversely,  the  Jones  vector  can  also  be  derived  from  the  polarization 
ellipse  using  the  general  equation 


E(x,y)  ~ 


■jS 

cos  <p 

-sin  q> 

COS  T 

sin  q> 

cos  <p 

j  sin  t 

(3-6) 


2.  Polarization  States 

The  Jones  vector  conveniently  describes  many  different  polarization 
states  -  especially  when  further  simplification  is  performed.  The  phase  8X  is  set 

to  zero  as  a  reference,  so  that  8y  =  8 .  Then  the  Jones  vector  can  be  expressed 
in  the  form  [12]: 

77  -  ^Ox 

0  ~  e  eJS 
_noye 

Using  Euler’s  formula  we  write 

be'5  =  b{  cos  8  +  j  sin  8)  =  B  +  iC  (3-8) 

so  that  the  Jones  vector  for  this  case  is 


A 

beiS 


(3-7) 


Eo  = 


A 

B  +  jC 


(3-9) 
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This  form  of  the  Jones  vector  is  completely  general  and  represents  all 
polarization  states  through  variations  of  the  parameters  A,  B  and  C.  A  normalized 

form  exists  whereby  the  Jones  vector  is  divided  by  V^42  +52  +C2  . 
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B.  PARTIAL  POLARIZATION  AND  STOKES  VECTOR 

So  far,  we  have  considered  only  monochromatic  waves  that  are 
completely  polarized  and  in  which  the  other  wave  attributes  (amplitude  and 
phase)  are  time  invariant.  A  discussion  of  partial  polarization  is  included  here  for 
completeness  and  to  illustrate  the  usefulness  of  Stokes  vectors  representations. 

By  partial  polarization,  we  mean  a  quasi-monochromatic  wave  with  time- 
dependent  amplitude  and  phase.  The  wave  polarization  traces  an  ellipse  with 
time-varying  tilt  and  amplitude.  As  further  illustration,  Mott  [11]  characterized  a 
quasi-monochromatic  wave  to  be  a  sum  of  a  completely  polarized  and  a 
completely  unpolarized  wave  from  which  the  degree  of  polarization  can  be 
defined. 

1.  Coherency  Matrix 

Using  this  decomposition  concept,  and  because  of  the  time-variant  or 
stochastic  nature  of  partially  polarized  waves,  a  coherency  matrix  containing  the 
correlations  between  polarized  and  unpolarized  components  of  the  electric  field 
is  used.  The  derivation  of  the  coherency  matrix  (as  detailed  in  [2]  and  [11])  can 
be  visualized  through  geometric  optics  because,  optically,  coherency  is  the  ability 
to  form  interference  fringes  [2].  Recalling  the  definition  of  the  E-field  from  (2.12), 
E  =  Re[E0ej{ox~h) ]  and  propagating  this  E-field  through  a  bire-fringent  phase  plate 

that  introduces  a  phase  difference  8  between  the  x  and  y  components,  we  may 
re-write  the  field  as  one  propagating  along  the  z-axis  and  making  an  angle  6  to 
the  x-axis,  so  that: 

E0  (t,  6)  =  Ex  cos  6  +  EyeJS  sin  0  (3.10) 
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Figure  111-1  Illustration  of  Wave  Propagating  in  z-direction 
The  time-averaged  intensity  is  then: 

<1^0  1 2>=Jxxcos20  +  Jyysm20  +  (Jxye-jS+Jyxejd)cos0sin0  (3.11) 

where  [J]  is  the  coherency  matrix  and  the  elements  are  averaged  over  a  number 
of  oscillations  or  cycles  as  denoted  by  the  angular  brackets: 
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Unpolarized  waves  have  time  averages  that  are  independent  of  the  angle 
the  wave  makes  to  the  x-y  axis  as  it  propagates  along  the  z-direction;  i.e.,  over 
time,  the  amplitude  of  the  power  or  time-averaged  intensity  does  not  depend  on 
the  polarization  of  the  E-field  since  the  E-field  has  no  constant  polarization  state. 
This  requires  that: 

J„=J„  (3.13) 

Jv=0  =  J„  (3.14) 

The  coherency  matrix  is  important  in  partial  polarization:  its  trace, 
determinant  and  eigenvalues  all  have  a  direct  mathematical  relationship  to  the 
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degree  of  polarization,  P.  The  coherency  matrix  can  also  be  subjected  to  unitary 
transform  so  that  it  can  be  decomposed  into  its  unpolarized  and  polarized 
components.  However,  the  coherency  matrix  is  difficult  to  visualize  and, 
alternatively,  Stokes  vectors  can  be  used  to  characterize  the  amplitude  and 
polarization  of  partially  polarized  waves. 


2.  Stokes  Vector 

For  completely  polarized  or  monochromatic  waves,  the  Stokes  vectors  are 
introduced  simply  as  [11]: 


S0  =  \Ex\2+\Ey\2 

(3.15) 

CN _ 

1 

<N 

II 

(3.16) 

S2=2\Ex\\Ey\cos8 

(3.17) 

S3=2\Ex\\Ey\sm8 

(3.18) 

where  E  and  8  are  as  previously  defined.  It  can  be  seen  that  from  the  Stokes 
vectors,  all  information  about  the  amplitude  and  polarization  can  be  defined.  In 
particular,  note  also  that: 

Sl=Sl+Sl+S\  (3.19) 


The  Stokes  vector  is  often  expressed  in  matrix  form  as: 


S3 


(3.20) 


We  can  also  relate  the  Stokes  vectors  to  the  polarization  ellipse  parameters: 


tan  2q>  = 


cos  8  =  — 


(3.21) 
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sin  2r  =  2'E*'Ey  sinS  =  ^- 

\EX\  2  +  Ey  S0 


(3.22) 


In  the  partially  polarized  case,  we  separate  the  Stokes  vectors  into  the 
polarized  and  unpolarized  components  as  before, 


(3.23) 
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With  some  manipulation  we  establish  the  relationship  between  the  coherency 
matrix  and  the  Stokes  vector. 


Sr,  =  J  +J 

0  xx  yy 

(3.24) 

e  —  /  _  t 

°1  U  XX  u  yy 

(3.25) 

e  —  /  i  / 

J2  u  xy  ^  u  yx 

(3.26) 
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II 

(3.27) 

Then  for  the  completely  unpolarized  case,  from  (3.13)  and  (3.14),  we  have 


Sl  =  S2  =  S3  =  0 


and, 


(3.28) 


~RS0~ 

"(l  ~R)S0~ 

Si 

0 

S  = 

+ 

0 

_ 

0 

(3.29) 


where  R  is  the  degree  of  polarization  and  can  be  defined  as: 


R  = 


S0pQl  yj +  S’j  +  51; 


=  1- 


4det{[J]} 


(3.30) 
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Note  that  for  the  completely  polarized  case,  R=1  and  the  equation  reduces 
back  to  (3.20). 


3.  Poincare  Sphere 

The  Poincare  sphere  is  a  geometrical  visualization  of  the  Stokes  vectors. 
Re-expressing  the  Stokes  vectors  in  terms  of  ^ ,  <p  and  r  (from  the  polarization 
ellipse),  we  obtain  [1 1]: 


51  =  S0  cos  2 t  cos  2<p 

52  =  S0  cos  2 t  sin  2 (p 


(3.31) 

(3.32) 


S3  =  S0  sin  2 t 


(3.33) 


It  is  then  apparent  that^ ,  S2  and  S3  are  the  Cartesian  coordinates  of  a  point  on  a 
sphere  with  radius  S0.  2cp  and  2r  are  the  azimuthal  and  elevation  angles 

respectively.  It  is  also  popular  to  express  the  Stokes  vectors  in  terms  of  the 
Deschamps  parameters  that  are  mathematically  defined  as  [8]: 


'Sy' 

r  cos  2  y 

S2 

=  So 

sin2/cos2£- 

vsin2^sin2£'y 

(3.34) 


Hence,  any  polarization  state  can  be  described  by  a  point  on  the  sphere  and  vice 
versa. 
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In  addition  to  the  Stoke’s  vectors,  the  Poincare  sphere  also  carries 
information  about  the  p,  q  and  P  (the  various  polarization  ratios  defined  in 
Chapter  II)  that  can  be  visualized  on  the  sphere  itself.  This  is  done  by  mapping 
the  Poincare  sphere  onto  the  complex  plane.  Before  doing  so,  we  summarize  the 
parameters  for  some  cases  of  polarizations. 

Table  111-2  Polarization  Parameters  for  Some  Polarization  States  [11] 


Polarization 

P 

p 

q 

$1 

s2 

S3 

Right  Circular 

1 

-j 

0 

0 

0 

-So 

Left  Circular 

-1 

j 

00 

0 

0 

So 

Linear  Vertical 

j  00 

00 

-1 

So 

0 

0 

Linear  Horizontal 

0 

0 

1 

1 

So 

0 

Figure  111-3  shows  the  way  to  map  a  point  on  the  Poincare  sphere  onto  the 
p,  P  complex  planes  to  obtain  their  values.  The  north-pole  (z-axis)  of  the  sphere 
is  aligned  with  the  real  axis  of  the  p  plane  and  the  equator  of  the  sphere  is 
aligned  with  the  imaginary  axis.  The  mapping  is  done  by  drawing  a  line  that 
intersects  (i)  the  furthest  point  on  the  sphere  from  the  plane  (specifically  the 
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highest  point  along  the  x-axis  in  this  illustration),  (ii)  the  point  of  interest  in  the 
sphere  and  (iii)  the  projection  point  on  the  p  complex  plane. 

In  Figure  111-3,  point  A,  which  is  also  the  north  pole  of  the  sphere, 
describes  p=-1  and  P=  j,  which  is  a  projection  from  the  point  (0,  0,  S0)  of  the 

sphere  (with  Cartesian  coordinates  (S', ,  S2 ,  S3)).  Point  B  describes  p=j  and  P=1, 
which  is  projected  from  the  point  (0,  S0 ,  0)  from  the  sphere.  This  construction 

agrees  with  Table  111-2.  Any  other  point  on  the  sphere  referring  to  any  polarization 
states  can  be  mapped  in  a  similar  manner. 


Figure  111-3  Mapping  of  Poincare  sphere  onto  p  and  P  Complex  Plane  [11] 


C.  SCATTERING/  SINCLAIR  MATRIX 

When  an  EM  wave  is  incident  on  a  target,  the  polarization  of  the  back- 
scattered  wave  is  generally  different  from  the  incident  wave  and  the  geometry  of 
the  target  will  determine  the  extent  of  depolarization.  Sinclair  proved  that  the 


30 


target  acts  as  a  polarization  transformer  and  defined  this  change  by  a  2  x  2 
scattering  matrix  [8]. 


[A]  = 


(3.35) 


The  important  concept  of  the  scattering  matrix  is  that  it  links  the  two 
orthogonal  components  of  the  incident  Jones  vector  to  the  back-scattered 
components. 


~K 

l 

A 

A 

_K_ 

4Arr 

A 

A 

Ey_ 

(3.36) 


The  distance  term  _  is  often  omitted  from  literature  since  only  the 

V4  nr 

polarization  is  of  interest  most  of  the  time.  For  a  mono-static  system  (where  the 
radar  transmitter  and  receiver  are  co-located),  the  reciprocity  theorem1  states 
that  the  cross  polarization  terms  are  equal  [9]. 

Ay  =  4.  (3-37) 


1.  Examples  of  Scattering  Matrices  of  Canonical  Targets 

The  scattering  matrices  are  dependent  on  the  geometry  of  the  targets. 
Once  the  scattering  matrix  is  specified,  the  back  scattered  wave  can  be 
computed  for  any  incident  wave  polarization  state.  The  scattering  matrices  of  a 
number  of  canonical  targets  have  been  well  established  in  literature  and  are 
tabulated  here. 


1  Reciprocity  Theorem  -  If  a  voltage  source  E  acting  in  one  branch  of  a  network 
causes  a  current  I  to  flow  in  another  branch  of  the  network,  then  the  same 
voltage  source  E  acting  in  the  second  branch  would  cause  an  identical  current  I 
to  flow  in  the  first  branch. 
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Table  111-3  Scattering  Matrices  of  some  Canonical  Targets  [8], [2] 


2.  Determination  and  Measurement  of  the  Scattering  Matrix 

The  scattering  matrix  contains  all  the  information  about  the  scattering 
properties  of  the  target.  The  elements  of  the  scattering  matrix  are  generally 
complex  quantities  that  may  be  expressed  as  [14] 


[A]  = 


I A  1 

\eiK 
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\e*"~ 

|  XX  | 

\  %y  \ 
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\e*» 
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_\  yx\ 

yy\ 

1 

(3.38) 


The  components  of  the  matrix  are  related  to  the  scattering  cross  section 
by 

<T*=4^2K|2>  axy  =  4^rr2 
Gyx=4nr2\Ay^ ,  <7^=4  nr1 

where  r  is  the  range  of  the  target.  Hence  the  specification  of  the  scattering  matrix 
is  dependent  on  the  measurement  of  the  scattering  cross  section  of  the  target.  In 
the  mono-static  case  A  =  Ayx ,  and  this  leaves  three  amplitudes  and  three 

phases  to  evaluate.  The  absolute  phase  depends  on  the  distance  of  the  target  to 
the  radar  and  the  target  orientation.  It  is  almost  impossible  to  measure  the 


xy 


yy 


(3.39) 
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absolute  phase  accurately  and  a  relative  phase  measurement  is  used  instead.  In 
this  case,  one  of  the  phase  terms,  e.g.,  ^xx,  is  factored  out  and  the  scattering 
matrix  becomes 


ere 


l(Tvye 


■sf\xr 


(3.40) 


Then,  five  quantities  are  left  to  be  measured  for  the  monostatic  case:  three 
amplitudes  and  two  relative  phases.  Relative  phases  are  easier  to  measure  than 
the  absolute  phase  since  they  are  not  dependent  on  the  target  distance.  Also, 
absolute  phase  is  not  necessary  for  most  cases  to  determine  the  polarization 
states.  Ruck  et.  al.  [14]  provides  the  elements  to  be  measured  in  a  laboratory 
transmit  and  receive  scenario  for  both  absolute  and  relative  phase  method. 


3.  Optimal  Polarizations 

We  can  make  use  of  the  scattering  matrices  and  the  polarization  match 
factor  or  polarization  ratios  introduced  in  Chapter  II  to  determine  the  polarization 
whereby  the  received  power  is  a  maximum  or  a  minimum. 

The  scattering  matrix  in  one  polarization  basis  may  be  transformed  to 
another  polarization  basis  in  a  different  orientation,  e.g.,  in  a  bi-static  scenario 
whereby  the  transmitting  and  the  receiving  antenna  may  not  have  the  same 
polarization  states.  This  is  done  using  a  unitary  transformation  matrix  [4], [8]: 


where  , 


cos  a 
-  sin  a  e 


-is 


sin  a  ejS 
cos  a 


(3.41) 


(3.42) 


and  a  describes  the  change  in  orientation  angle  between  the  two  polarization 
bases.  We  have  previously  defined  the  polarization  ratio  in  (11-21)  as  P  =  tan<p  eJS 
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relating  the  orthogonal  components  of  an  EM  wave.  It  can  be  redefined 
analogously  as  a  polarization  ratio  that  describes  the  change  in  polarization  basis 
as: 


P  =  tan  a  ejS 


(3.43) 


In  this  case,  the  unitary  transformation  matrix  can  then  be  written  simply  as: 


P] 


1 

Jl  +  PP* 


P 

1 


Now  the  elements  of  the  [A\b  matrix  can  be  written  as: 


Aa  =  (1  +  PP'rl[A„  +  P2AW  +  P(AV  +  AJ] 

Aab  =  (1  +  PP'  T\-P'A„  +  PAW  +Axy-  PP*  A  yx)] 
Aba=d  +  PPri[-P*Axx+PA>y+Ayx -PP*Axy)\ 
Abb  =  (1  +  PP*)~l  [-P*2Axx  +Ayy-P*  (Axy  +  Ayx )] 


(3.44) 


(3.45) 


In  the  optimization  problem,  we  can  visualize  the  two  scattering  matrices 
as  the  transformation  in  polarization  state  when  the  EM  wave  is  depolarized  by  a 
target.  The  power  of  the  scattered  wave  is  distributed  into  two  channels  of 
reception,  the  co-polarized  channel  (identical  to  the  polarization  state  of  the 
transmitting  antenna)  and  the  cross-polarized  channel  (orthogonal  to  the 
polarization  state  of  the  receiving  antenna). 


Figure  MI-4  Polarization  States  of  Transmitted  and  Backscattered  Waves 

It  was  first  shown  by  Kennaugh  that  there  exist  two  pairs  of  optimal 
polarizations  states  [1].  The  co-polarization  null  pair  (CO-POL)  for  minimal  co¬ 
polarization  (and  maximum  cross-polarization)  is  obtained  by  setting  Aaa  and  Abb 
to  zero  so  that  [2]: 
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(3.46) 


The  cross  polarization  null  pair  (X-POL)  for  maximal  polarization  is  obtained  by 
setting  Aab  and  Aba  to  zero  so  that  [2]: 


px  _  ~R\  ±  A2  +  4 ^2 ^3 

1,2  2  R2 

=(\Aj-\ 4,  |2  -AxyA*yx  +  AyxA*J  (3.47) 

*2=(VC  +  VO 

R3=(AXXAly+AxyAly) 

The  CO-POL  and  the  X-POL  pairs  are  often  plotted  onto  the  Poincare  sphere 
forming  what  is  known  as  the  Polarization  fork,  originally  conceived  by  Huynen. 

Many  optimization  techniques  have  been  developed  and  a  good  summary 
is  found  in  [1],  Such  techniques  may  involve  matrix  manipulation  of  a  target 
scattering  matrix  to  derive  its  eigenvalues,  deriving  the  polarization  ratios  from 
the  Poincare  sphere  or  using  the  Stokes  vector  and  Muller  matrix.  Optimization  of 
polarization  states  will  be  further  discussed  in  Chapter  IV  through  its 
implementation  in  Matlab. 


D.  MUELLER  (POWER)  MATRIX 

For  the  partially  polarized  case,  it  is  necessary  to  relate  the  incident  and 
scattered  Stokes  vector.  Just  as  the  incident  and  scattered  Jones  vector  is 
related  by  the  scattering  matrix,  the  incident  and  scattered  Stokes  vector  is 
related  by  the  Mueller  matrix  [8]: 


J  scattered 


=  [M]Su 


incident 


(3.48) 


We  express  the  Stokes  vector  in  the  energy  form  to  give  [9]: 
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(3.49) 
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Then  relation  between  the  incident  and  scattered  energy  is  given  by  the  relation 


(3.50) 


so  that  [14], 
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(3.51) 

Generally,  the  scattering  matrix  contains  complex  elements  of  a  voltage  and  is 
used  for  coherent  systems  while  the  Mueller  matrix  contains  power  elements  and 
is  used  for  incoherent  systems[14]. 


E.  TARGET  DECOMPOSITION  THEOREMS 

Cloude  [5]  defines  the  main  idea  behind  Target  Decomposition  (TD) 
theorems  as  “expressing  the  average  scattering  matrix  for  a  random  media 
problem  as  a  sum  of  independent  elements  and  to  associate  a  physical 
mechanism  with  each  component.” 
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The  scattering  matrix  is  a  measured  quantity,  often  stochastic  rather  than 
deterministic  in  nature.  The  use  of  TD  theorems  to  evaluate  the  average 
scattering  matrix  allows  the  various  scattering  mechanisms  and  target  physical 
properties  to  be  identified.  In  addition,  it  is  possible  to  generate  a  dominant 
scattering  mechanism  for  better  physical  interpretation.  This  is  the  main 
advantage  of  TD  theorems. 

The  method  described  here  is  based  on  developments  by  Cloude  and 
involves  the  use  of  Pauli  matrices  to  vectorize  the  scattering  matrix.  The 
scattering  vector  can  then  be  related  to  the  physical  attributes  of  the  scattering 
target  using  the  Scattering  Vector  Reduction  Theorem.  Finally,  these  physical 
attributes  can  be  used  to  identify  the  target  against  a  cluttered  background  that 
has  different  scattering  mechanisms. 

1.  Pauli  Spin  Matrices 

The  Pauli  spin  matrices  (commonly  used  in  quantum  mechanics)  can  be 
used  to  decompose  the  scattering  matrices.  In  quantum  mechanics,  Pauli 
matrices  are  a  basis  set  that  described  microscopic  particles  of  spin  V2.  Spin  is  a 
quantity  similar  to  angular  momentum  in  classical  mechanics  but  is  distinguished 
by  its  quantized  and  discrete  nature  [19].  The  Pauli  spin  matrices  are  [9]: 

(3.52) 


The  interpretation  of  the  Pauli  matrices  is  as  tabulated: 
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Table  111-4:  Pauli  Matrices  and  their  Interpretation  [9] 


Pauli  matrix 

Scattering  Type 

Interpretation 

"1  0" 
0  1_ 

Odd-  bounce 

Surface,  sphere,  corner  reflectors 

"1  0 

0  -1 

Even-bounce 

Dihedral 

'0  f 
t  0_ 

Even-bounce  n/4  tilted 

n/4  tilted  dihedral,  volume  scattering 

'0  -j 

J  0 

Cross-polariser 

Non-existent  for  backscattering 

The  advantage  of  using  a  Pauli  matrix  basis  is  that  the  elements  that 
result  are  closely  related  to  the  physics  of  wave  scattering  [4].  The  scattering 
types  that  result  after  decomposition  are  orthogonal  so  that  their  separation  is 
straightforward.  As  an  example,  recall  the  concept  presented  in  Table  Table  111-3 
that  the  polarization  states  of  backscattering  from  targets  are  determined  from 
their  scattering  matrices.  We  illustrate  this  again  in  the  following  Figure  111-5 
(notice  that  the  two  scattering  matrices  form  two  of  the  Pauli  matrices  defined 
above). 


Figure  MI-5  Relating  Scattering  Matrices  to  Pauli  Matrices 
Based  on  this  idea,  the  Pauli  decomposition  can  be  used  to  process  a  full 
polarimetric  data  set  and  separate  them  into  their  orthogonal  sets.  Figure  MI-6  of 
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the  Great  Aletsch  glacier,  Switzerland  gives  a  vivid  illustration  of  this  property  of 
the  Pauli  matrices  [18].  Using  Pauli  decomposition,  a  three  color  picture  of  a  fully 
polarimetric  data  set  is  constructed.  Areas  with  even-bounce  characteristics 
represent  buildings  in  urban  areas  and  are  coded  in  red.  Areas  with  odd-bounce 
characteristics  representing  flat  surface  and  corner  reflectors  are  coded  in  blue. 
Lastly  areas  with  scattering  mechanism  of  the  third  Pauli  matrix  that  represents 
volume  scattering  are  coded  in  green. 


Figure  111-6  Pauli  Decomposition  of  a  Full  Polarimetric  Data  Set  [9] 


2.  Vectorization  of  Scattering  Matrices 

The  extraction  of  physical  information  about  a  target  involves  the 
vectorization  of  its  scattering  matrix.  The  vectorization  process  makes  use  of  the 
Pauli  matrices  as  basis  elements  to  separate  the  scattering  matrix  into  its 
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physical  scattering  mechanism.  The  representation  of  the  scattering  matrix  in  the 
Pauli  basis  matrix  results  in  a  vector  k  through  the  following  operation  [5], [7]: 
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In  the  case  of  backscattering,  the  reciprocity  theorem,  as  mentioned  in 
equation(3.37),  results  in  AHV=AVH,  thereby  eliminating  the  fourth  term  in  the 
scattering  vector,  reducing  it  to: 
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V2 
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^HH  ^  ^VV 
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^VV  ^HH 

2  Ahv 


(3.54) 


Specifically,  the  elements  of  this  scattering  vector  are  derived  from  a 
scattering  matrix  decomposed  by  a  Pauli  basis  set.  The  first  element  AHH  +  AVV  is 
from  the  first  Pauli  matrix  and  hence  represent  an  odd-bounce  scattering  type. 
The  second  element  Ayy-AHH  corresponds  to  an  even  bounce  scattering  type 

and  the  third  element  2 AHV  corresponds  to  the  third  Pauli  matrix.  This  then 

illustrates  how  the  Pauli  decomposition  of  a  full  polarimetric  data  set  is 
accomplished.  Referring  back  to  Figure  111-6,  areas  coded  blue  are  the  first 
elements  of  the  scattering  vector,  areas  coded  red  are  the  second  elements  and 
areas  coded  green  are  the  third  elements. 

In  order  to  correlate  the  target  scattering  vector  to  the  actual  physical 
attribute  and  scattering  mechanism,  we  need  to  equate  the  scattering  vector  to  a 
general  unitary  complex  vector  [7]  whose  parameters  are  linked  to  physical 
space,  so  that  now, 
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(3.55) 
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Equation  (3.55)  may  be  visualized  as  a  unitary  transformation  of  the  matrix 
to  a  different  axes  set  (e.g.  in  the  case  of  (3.36)).  One  way  to  do  this  is  to  utilize 

the  Poincare  sphere.  Hence  the  expression  of  the  k  vector  as  a  unitary  complex 
vector  is  analogously  equivalent  to  rotation  of  the  Poincare  sphere  where  a  and 
P  are  the  Deschamps  parameters,  previously  described  in  (3.34)  and  repeated 
here  [7]: 
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To  change  the  scattering  vector  k  into  k'  with  a  differential  change  in 
scattering  mechanism,  we  can  introduce  changes  Aa  and  Ap  in  the  angles  a  and 
p  respectively  [9].  Mathematically  this  is  done  by  using  vector  rotation  algebra: 
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This  observation  leads  to  the  important  concept  of  the  Scattering  Vector 
Reduction  Theorem  [4]  that  states  that  “it  is  always  possible  to  reduce  an 

arbitrary  scattering  mechanism  represented  by  a  complex  unitary  vector,  kto  the 
identity  [l,0,of  by  the  following  set  of  ordered  matrix  transformations”: 
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We  are  now  finally  ready  to  interpret  the  above  in  terms  of  both  the 
physical  attributes  of  the  target  as  well  as  its  scattering  mechanism.  The  two 
matrices  above  resemble  the  matrices  for  plane  rotations.  Physically,  p 
corresponds  to  the  rotation  of  the  sensor  coordinates  [4].  Conversely,  this  angle 
also  refers  to  the  orientation  of  the  target  and  takes  on  a  value  of  -180°  to  180°. 

The  first  matrix  containing  the  angle  a,  though,  also  resembles  a  plane 
rotation  and  is  not  related  to  target  orientation.  This  matrix  corresponds  to  an 
internal  degree  of  freedom  of  the  target  and  hence  is  directly  related  to  the  type 
of  scattering  mechanism  of  the  target  [4],  The  angle  a  is  continuous  within  a 
range  of  0°  to  90°  and  can  be  used  to  represent  a  variety  of  different  scatterers. 


anisotropic  odd  bounce  anisotropic  even  bounce 
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Figure  111-7  Interpretation  of  the  a  angle  [9] 


The  interpretation  of  the  a  angle  is  proposed  and  developed  by  Cloude  [6]. 
He  suggests  that  the  scattering  by  a  media  can  be  modeled  by  a  single  dominant 
matrix,  the  coherency  matrix,  defined  as 

[T]  =  k-k*T  (3.60) 

The  coherency  matrix  [T]  has  three  eigenvalues,  X,  and  and  the 

entropy,  H,  of  these  eigenvalues  is  a  measure  for  the  randomness  of  the 
averaged  polarimetry  scattering  with  H  =  0  indicating  a  single  scattering 
mechanism  and  H  =  1  representing  a  random  mixture  of  scattering  mechanisms, 


42 


i.e.  a  depolarizing  target.  Values  between  0  and  1  indicate  the  degree  of 
dominance  of  one  particular  scatterer: 
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The  angle  a  can  be  determined  by  the  parameterization  of  [r] 
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(3.64) 


a  =  Pxax  +  P2a2  +  P3a3 


(3.65) 


An  cc-H  plot  can  be  constructed  that  defines  9  regions  for  different  scattering 
zones  at  different  entropies  (ranging  from  Low  Entropy  Surface  Scattering  to 
High  Entropy  Multiple  Scattering). 
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Figure  111-8  Zones  in  a-H  plane  for  Random  Media  Scattering  Problems  [6] 

At  a=  0°,  we  have  an  isotropic  surface  represented  by  the  first  Pauli 
matrix.  As  a  increases,  the  surface  becomes  increasingly  anisotropic.  At  a=90°, 
we  obtain  the  second  Pauli  matrix.  This  point  can  also  be  used  to  represent 
targets  that  cause  a  phase  shift  of  2n  between  the  H  and  V  components  of  the 
electric  field.  For  intermediate  values  of  a ,  the  scattering  mechanism  becomes 
anisotropic  and  AHH  and  A  are  no  longer  equal.  The  anisotropic  behaviour  and 

consequently  the  increasingly  different  co-polar  scattering  coefficients  means  a 
preferential  scattering  in  a  particular  direction.  At  the  limit  of  a=45°,  this  leads  to  a 
dipole  scatterer  whose  orientation  is  determined  by  the  angle  p  and  one  of  the 
copolarized  coefficient  goes  to  zero. 

The  a  angle  is  decoupled  from  the  p  angle  so  that  scattering  mechanism 
can  be  identified  independently  of  the  target’s  physical  orientation.  A  list  of  the 
possible  a  and  p  angles  are  tabulated  and  mathematical  examples  of  their 
derivation  are  given  in  Chapter  IV. 
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Table  MI-5  Examples  of  a  and  p  angles  for  Canonical  Targets  [9] 


Canonical  Scatterer 

a 

p 

Sphere 

0° 

-180°  to  180° 

Dihedral  at  0 

CD 

O 

o 

20 

Dipole  at  0 

45° 

20 

Helix 

CD 

O 

o 

±45° 

Surface  at  0 

0° 

20 

Many  other  TD  theorems  exist  and  different  approaches  seek  to  extract 
different  information.  The  TD  theorems  may  also  be  applied  to  the  various 
polarization  matrices  besides  the  scattering  matrix  such  as  the  coherence  and 
the  Mueller  matrices.  A  good  review  of  the  various  approaches  is  available  in  [5] 
and  [9],  pg  43-58. 
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IV.  POLARIMETRY  THEORETICAL  ANALYSIS 


The  approach  of  this  analysis  is  to  use  numerical  examples  to  verify  the 
TD  theorems  put  forth  by  current  researchers.  The  objective  is  to  illustrate  that 
full  polarimetric  data  set  does  enhance  target  detection,  simply  by  the  fact  that 
target  information  is  lost  in  the  target  returns  through  the  cross-pol  channels  that 
are  not  collected.  Depending  on  the  target  geometry  and  orientation,  the 
recovery  of  such  information  may  be  pivotal  to  the  detection  of  the  target. 

The  second  aspect  of  the  analysis  is  to  create  a  Matlab  model  that  can 
derive  the  optimal  polarization  states  of  the  incident  electric  field  for  any  target 
scattering  matrix.  There  currently  exist  many  analytical  techniques  to  obtain  the 
optimal  polarization  states  and  our  implementation  in  Matlab  provides  a  fast  and 
easy  numerical  instantiation  of  the  derived  solutions. 

The  source  codes  for  the  two  models  are  included  in  the  appendices. 


A.  USING  TD  THEOREMS  FOR  TARGET  DETECTION 


1.  Scattering  Matrices  and  Vectors 

Before  we  use  the  TD  theorems  in  Chapter  III,  some  verification  was 
conducted  to  ensure  commutability  between  the  scattering  matrix  and  scattering 
vector  through  known  matrices  of  canonical  targets,  also  introduced  in  Chapter  3. 

Expressing  equation  (3.52)  in  terms  of  the  elements  of  the  scattering 
matrix  yields 
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(4.1) 


Starting  with  the  sphere  or  the  flat  plate,  the  angle  a  is  selected  as  0°  for 
an  isotropic  odd  bounce  surface  while  J3  can  adopt  any  angle  due  to  the 
symmetry  of  the  sphere.  Substituting  this  into  the  equation  (4.1)  we  obtain 
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(4.2) 
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or,  in  the  more  common  scattering  (normalized)  form, 
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This  is  in  agreement  with  Table  111-3  of  scattering  matrices  of  canonical  targets. 

Similarly  for  a  dihedral  orientated  at  an  angle  6 ,  the  a  angle  is  selected 
as  a  90°  for  an  isotropic  even  bounce  surface,  the  scattering  matrix  elements  are, 


f 


v 


A  \ 
ahh 

Arv 

A 

^HV  J 


2 


^cos  90  +  sin  90  cos  2d 
cos  90  -  sin  90  cos  20 
,  sin  90  sin  2# 


_  1 

r  cos  20  N 

— cos  20 

~2 

,  sin  2  0  , 

(4.4) 
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cos  20  sin  20 
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(4.5) 


The  above  examples  illustrates  that  the  target  scattering  vector  does 
connect  the  scattering  matrix  to  the  target  geometry  and  scattering  mechanism 
by  appropriate  selection  of  the  a  and  p  angles. 


2.  Construction  of  Foliage 

Foliage  commonly  constitutes  radar  volume  clutter  and  hides  legitimate 
targets.  It  is  an  important  aspect  of  radar  detection  to  study  and  model  this 
target/clutter  interaction.  The  foliage  is  expected  to  bring  about  some  degree  of 
depolarization  of  the  electric  waves  but  this  effect  is  not  significant  in  any  one 
single  direction  due  to  the  random  distribution  and  orientation  of  the  leaves. 

We  use  the  same  basic  equation  derived  from  TD  analysis  and,  as  before, 
model  the  foliage  as  a  two-dimensional  layer  of  leaves.  The  model  is 
implemented  in  MatLab.  Each  leaf  resides  in  a  single  cell  in  a  grid  of  user-defined 
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dimensions  (We  use  a  50  by  50  grid  in  this  simulation).  Each  cell  is  then 
represented  by  a  scattering  vector  or  matrix.  This  concept  is  shown  in  Figure  IV-1 


cos  an 

y 


\ 


sincr^cos  /?y 
vsina..sin/?..y 


Figure  IV-1  Foliage  Modeling 

The  angle  a  is  selected  as  45°  to  simulate  the  leaves  as  a  cloud  of  small 
anisotropic  dielectric  dipoles.  The  p  angle  is  taken  to  be  uniformly  distributed 
from  0  to  n  to  simulate  the  random  orientation  of  the  leaves.  From  the  scattering 
vectors,  we  derive  the  scattering  behavior  and  populate  each  cell  with  an 
appropriate  scattering  matrix. 


Now,  we  introduce  a  propagating  electric-field  to  interact  with  the  foliage 
and  then  study  the  backscattered  polarization  states.  Jones  vectors  are  used  to 
represent  horizontally,  vertically  and  circularly  polarized  E-fields. 


Table  IV-1  Jones  Vectors  of  Incident  E-fields 


Horizontal  Linear 

Polarization 

£„  = 
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Linearly  Polarized  at 
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Eo  = 

"  l 

71 
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Left  Circular 

Polarization 
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°  M 

T 
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The  Jones  vector  will  interact  with  each  of  the  cells  according  to  the 
following  equation: 
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The  output  of  our  model  is  a  graphical  display  of  the  horizontal  and  the 
vertical  components  of  the  backscattered  E-field. 

a)  Horizontal  Incident  E-Field 

The  plot  shown  below  is  the  backscattered  E-field  due  to  a 
horizontally  polarized  incident  E-field,  E'  =  [l  0]r . 


Horizontal  Scattered  E-field 


Vertical  Scattered  E-field 


Figure  IV-2  Backscattered  E-field  due  to  a  Horizontally  Polarized  Incident  E-field 
The  mean  scattering  matrix  of  the  foliage  from  this  simulation  run 
is: 


A 


foliage 


0.3514  0.2274 
0.2274  0.3558 
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and  the  mean  return  of  the  electric-field  is 


Er 


0.3514 

0.2274 


These  results  show  that  some  degree  of  depolarization  has  occurred  as  the 
completely  horizontally  polarized  incident  e-field  is  now  depolarized  by  the  foliage 
with  backscattered  E-field  components  in  both  vertical  and  horizontal  polarized 
states. 


b)  Linearly  Polarized  at  45°  Incident  E-field 

The  simulation  run  was  repeated  with  linearly  polarized  at  45°  E- 
field  with  components  of  equal  magnitude  in  both  horizontal  and  vertical 
polarization  state. 


Horizontal  Scattered  E-field 


Vertical  Scattered  E-field 


Figure  IV-3  The  Backscattered  E-field  from  a  Linearly  Polarized  at  45°  E-field 
The  mean  scattering  matrix  is  identical  to  the  last  simulation  and 
the  mean  backscattered  E-field  is  now: 
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E 


0.41 

0.40 


This  result  is  to  be  expected  since  the  incident  E-field  is  equally  polarized  in  the 
vertical  and  horizontal  directions  -  the  foliage  being  a  random  depolarizer,  it 
should  not  have  significant  effect  on  the  E-field  polarization  over  sufficiently  large 
sample. 

c)  Circularly-Polarized  Incident  E-field 

Another  interesting  scenario  that  was  examined  used  a  circularly 


polarized  incident  electric  field  with  the  Jones  vector, 


E  = 


V2’-'V2 


for  a  left 


circularly  polarized  incident  field.  The  results  for  this  case  are  presented  in  Figure 
(V-4)  (both  the  real  and  imaginary  parts  of  the  echo  field). 


Figure  IV-4  Backscattered  E-field  from  a  Left  Circularly  Polarized  Incident  E-field 
The  mean  of  the  backscattered  E-field  is  then 
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(4.7) 


0.25+  /  0.16 
0.16  +  y'0.25 


3.  Modeling  Target  and  Foliage 

The  idea  now  is  to  locate  a  target  within  the  foliage  and  observe  the 
overall  scattering  behavior.  A  number  of  canonical  targets  may  be  used:  e.g.  flat 
plates,  dihedrals  oriented  at  different  angles  or  metallic  helixes.  The  scattering 
vectors  of  the  targets  are  listed  in  Table  MI-3.  The  location  of  the  target  is 
randomized  in  the  50  by  50  grid  so  that  each  simulation  run  is  associated  with  a 
different  target  location.  Figure  IV-5  gives  an  illustration  of  the  model. 


Dihedral  oriented  at 
an  angle  located  at  a 
random  location 


Figure  IV-5  Modeling  the  Target  and  the  Foliage  Together 


The  model  currently  holds  scattering  matrices  of  most  known  canonical 


targets  (including  linear  polarization  at  various  orientations,  dihedrals  at  various 


orientations,  flat  plates/  spheres  and  helixes)  and  Jones  vectors  of  various 


polarizations  of  the  incident  E-fields. 


The  backscattered  returns  will  now  be  due  to  both  the  foliage  clutter  and 
the  target  based  on  the  following  relationship 


Er  =  Am%etE  +  AclutterE  (4.8) 

a)  Horizontally  Polarized  Incident  E-field  on  Flat  Plate 

A  relatively  straight  forward  simulation  run  was  conducted  first 
based  on  the  scenario  of  a  horizontally  polarized  electric  field  impinging  onto  a 
flat  plate  target.  The  flat  plate  does  not  depolarized  the  incident  electric  field  and 
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the  reflected  electric-field  from  the  target  alone  is  therefore  still  expected  to  be 
located  in  the  horizontal  component.  The  result  of  this  simulation  run  is  given  in 
Figure  IV-6. 


Figure  IV-6  Backscattered  E-field  from  a  Horizontally  Polarized  Incident  E- 
field  reflected  from  a  Flat  Plate 


As  discussed,  for  this  case  the  bright  target  return  is  not 
depolarized  and  is  located  within  the  horizontal  components  of  the  backscattered 
E-field,  which  is  the  same  polarization  of  the  transmitting  antenna  that  produced 
the  horizontally  polarized  incident  field.  The  foliage  behaves  by  altering  the 
horizontal  incident  E-field  so  that  the  backscattered  E-field  has  a  vertical 
component.  However,  this  behavior  would  be  different  if  a  depolarizing  target  is 
used. 
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b)  Horizontally  Polarized  E-field  on  45 0  Dihedral  Target 

For  the  next  simulation  run  a  45°  oriented  dihedral,  known  to 
completely  depolarize  the  incident  field,  is  examined.  We  continue  to  use  a 
horizontally  polarized  incident  field.  The  result  obtained  is  shown  in  Figure  IV-7 


Horizontal  Scattered  E-field 


Figure  IV-7  Backscattered  E-field  from  a  Horizontally  Polarized  Incident  E-field 
reflected  from  a  45°  Oriented  Dihedral 

Now  the  target  return  is  completely  located  in  the  vertical 
component  of  the  backscattered  electric  field.  Hence  without  a  full  polarimetric 
data  set,  this  target  will  not  be  detected  despite  producing  such  a  bright  (clutter- 
free)  return. 

c)  Left-Circularly  Polarized  E-field  on  Flat  Plate  Target 

The  effect  of  circularly  polarized  incident  E-field  on  various  targets 
was  also  investigated  using  the  model.  The  Jones  vector  of  the  incident  E-field  is 
now 
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(4.9) 


1 


.  1 


and  it  represents  a  left-circularly  polarized  wave.  The  imaginary  vertical 
component  means  that  the  two  orthogonal  components  are  90°  out  of  phase.  The 
first  scenario  once  again  involves  a  flat  plate  target  and  the  result  of  the 
simulation  is  shown  in  Figure  IV-8  (including  both  the  real  and  the  imaginary 
components).  The  result  essentially  shows  that  the  backscattered  wave  is  not 
depolarized,  and  continues  to  be  a  circularly  polarized  wave  (with  real 
component  in  the  horizontal  and  imaginary  component  in  the  vertical). 


Horizontal  Scattered  E-field  Vertical  Scattered  E-field 


Imaginary  Horizontal  Scattered  E-field  Imaginary  Vertical  Scattered  E-field 


Figure  IV-8  Backscattered  E-field  from  a  Left-Circularly  Polarized  Incident  E-field 

reflected  from  a  Flat  Plate 


56 


d)  Left-Circularly  Polarized  E-field  on  Helix  Target 

The  last  scenario  included  here  is  that  of  a  circularly  polarized  E- 
field  impinging  on  a  helical  target.  The  helical  target  acts  as  a  polarizer  that 
brings  about  a  phase  difference  of  90°  to  the  vertical  and  horizontal  components 
of  the  incident  E-field.  Because  the  incident  E-field  components  are  already  out 
of  phase  by  90°,  the  helix  target  actually  removes  the  phase  difference  resulting 
in  a  linearly  polarized  backscattered  signal.  This  effect  is  shown  in  Figure  IV-9. 
The  backscattered  signal  is  now  a  horizontally  polarized  wave. 


Figure  IV-9  Backscattered  E-field  from  a  Left-Circularly  Polarized  Incident  E-field 

reflected  from  a  Helix 

This  portion  of  the  analysis  shows  the  possibility  of  using  TD 
theorems  to  extract  physical  properties  of  targets  from  their  scattering  matrices 
and  these  physical  attributes  can  be  used  in  modeling  -  in  this  particular  case, 
for  simulation  of  foliage.  The  latter  example  of  a  man-made  target  buried  in 
foliage  clutter  demonstrates  the  importance  of  obtaining  full  polarimetric  data  set 
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for  target  detection.  We  chose  a  particularly  bright  target  and  showed  that,  due  to 
depolarization,  the  return  signal  can  be  completely  lost.  It  should  also  be  noted 
that  for  the  dihedral,  the  characteristic  to  depolarize  also  depends  on  its 
geometry,  i.e.,  its  orientation  angle  of  45°. 

B.  POLARIZATION  OPTIMIZATION 

The  theory  for  optimization  of  polarization  state  is  detailed  in  [1]  and  a 
number  of  techniques  are  discussed  there.  For  simulation  in  Matlab,  we  select  an 
optimization  technique  that  uses  the  Stokes  vectors  and  the  Mueller  (power) 
matrix.  This  technique  has  the  particular  advantage  that  the  input  arguments  are 
in  the  power  or  intensity  units,  so  that  the  output,  which  is  also  in  terms  of  power 
units,  can  be  directly  obtained  by  selecting  the  appropriate  elements  of  the 
Mueller  matrix.  Since  the  program  will  essentially  vary  the  input  parameters  until 
the  peak  power  return  is  obtained,  the  intermediate  mathematical  techniques 
covered  in  [1]  can  be  avoided  entirely. 


1.  Modeling  the  Varying  Incident  Electric  Field  and  Stokes 
Vectors 

The  polarization  states  of  the  incident  electric  field  can  be  described 
complete  by  its  tilt  and  ellipticity  angles  (refer  to  Chapter  2  Section  B).  Equation 
(3.6)  shows  the  relationship  as  [8] 


cos©  -sin©  cos t 

E, .  =  Ae~jS 


sin  cp  cos  (p  II  j  sin  r 


(4.10) 


Flence,  any  Jones  vector  of  the  electric  field  may  be  described  by  varying 
the  tilt  and  the  ellipticity.  Flowever,  as  defined  by  the  polarization  ellipse,  the 
ellipticity  has  a  range  of  -45°  to  45°  while  the  tile  of  the  ellipse  has  a  range  of  0° 
to  180°  -  and  these  will  be  the  ranges  of  our  model  parameters. 

The  electric  field  is  used  to  generate  the  Stokes  vectors  using  the 
following  equation  [11] 
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(4.11) 


.  ,2  I  |2 

Ik  +KI 

l£  I2  -Le  I2 

I  x\  \  y\ 

2|£,|k|cos<J 


2\EA  Ey  sinS 


where  5  is  the  phase  difference  between  the  x  (horizontal)  and  y  (vertical) 
components  of  the  electric  field  and  can  be  derived  from  the  tilt  and/or  ellipticity 
angles  [12] 

\e  I2  —  \e  I2 

cos£  =  tan2<p — - -  —  (4.12) 

2|  Ex\ 


2.  Application  of  the  Mueller  (Power)  Matrix 

In  Chapter  III  we  explained  how  the  Mueller  matrix  transforms  the  incident 
Stokes  vectors  to  the  backscattered  Stokes  vectors,  analogously  to  the  operation 
of  scattering  matrix  and  Jones  vectors.  However,  we  need  to  separate  the 
backscattered  field  into  the  co-polarized  and  the  cross-polarized  channels  to 
explore  the  effects  of  polarization,  i.e.,  we  need  to  determine  how  much  of  the 
incident  power  is  depolarized  and  returns  in  the  orthogonal  component  (which  is 
collected  only  a  full  polarimetric  data  set). 

This  is  done  by  generating  two  Mueller  matrices,  one  each  for  the  co-pol 
and  the  cross-pol  channel,  in  a  manner  based  on  developments  by  Boerner  [1], 


M„  = 


10  0  0 
0  10  0 

0  0  10 

0  0  0  -1 


M 


M„  = 


10  0  0 
0-100 
0  0-10 
0  0  0  1 


M 


(4.13) 


(4.14) 
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Then,  the  co-pol  and  cross-pol  power  can  be  obtained  from  the  following 
equations  [1] 

p.MM!]  (4  15) 

INI  ‘  ’ 

(4  16) 

INI 


3.  Results  and  Analysis 

The  model  can  now  be  used  to  determine  the  polarization  state  for 
maximum  power  return  for  any  target  scattering  matrix.  The  choice  of  polarization 
state  of  the  incident  electric  field  ensures  a  maximum  return  in  either  the  co¬ 
polarized  or  the  cross-polarized  channel.  In  this  sense,  if  the  target 
characteristics  can  be  determined,  and  a  full  polarimetric  data  set  is  available, 
the  transmitted  electric  field  can  be  selected  to  optimize  the  power  return. 

The  first  scattering  matrix  used  for  simulation  was  that  of  a  flat  plate  (or  a 
sphere)  that  does  not  depolarize  the  incident  wave  field.  The  result  of  the 
simulation  is  presented  in  Figure  IV-10. 
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Figure  IV-10  Backscattered  Power  for  a  Flat  Plate 


The  figures  in  the  first  row  show  the  power  return  in  the  co-polarized 
channel  that  has  the  same  polarization  state  as  the  antenna  transmitting  and 
receiving  the  signal  (for  the  mono-static  case).  The  power  is  plotted  as  a  function 
of  the  tilt  and  ellipticity  angles.  The  figure  on  the  left  gives  a  three-dimensional 
presentation  while  the  right  figure  gives  a  two-dimensional  projection  in  which  the 
gradient  of  the  color  gives  the  power  return.  The  figures  in  the  second  row  give 
the  results  for  the  cross-polarized  channel. 


The  results  show  that  the  return  power  depends  only  on  the  ellipticity 
angle  and  not  the  tilt  angle.  For  a  0°  ellipticity  angle,  the  electric  field  is  linearly 
polarized  and  the  maxima  occurs  in  the  co-polarized  channel  showing  maximum 
power  is  backscattered  with  the  same  polarization  as  the  transmitting  antenna. 
As  the  ellipticity  increases  to,  say  10°,  the  co-pol  power  begins  to  drop,  similar  to 
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the  co-pol  power  of  the  incident  E-field.  When  the  ellipticity  is  increased  to  its 
maximum  of  45°,  the  cross-pol  power  attains  its  maximum  value. 


Figure  IV-1 1  Variation  of  Co-pol  and  Cross-pol  Powers  with  Ellipticity 
This  behavior  shows  that  the  target  does  not  depolarize  the  incident  field 
at  all.  This  result  is  expected  for  a  flat  plate  target  and  hence,  the  model  has 
accurately  simulated  this  behavior. 

The  next  scattering  matrix  used  for  simulation  was  that  of  the  45°  oriented 
dihedral  that  completely  depolarizes  the  incoming  electric  field.  The  appropriate 
scattering  matrix  is 

^45" dihedral  ~  1  f>  (4.17) 
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The  results  of  the  simulation  is  presented  in  Figure  IV-12  below. 


Figure  IV-12  Backscattered  Power  for  Dihedral  Oriented  at  45° 


The  results  given  in  Figure  IV-12  show  that  the  co-pol  nulls  are  the  cross- 
pol  peaks.  We  identify  the  tilts  and  ellipticity  whereby  the  cross  pol  peaks  occur. 

Table  IV-2  Incident  Electric  Fields  for  Cross-Pol  Peaks 


Tilt,  q> 

Ellipticity,  r 

Equivalent  Jones 

0° 

0° 

"0" 

1 

CO 

o 

o 

0° 

"f 

_0_ 

We  can  conclude  from  the  simulation  results  that  horizontally  and 
vertically  polarized  electric  fields  are  completely  depolarized  by  a  45°  oriented 
dihedral. 
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The  last  scattering  matrix  that  we  present  here  is 


A  = 


2  j 
0.5 


0.5 

~j 


(4.18) 


This  is  an  arbitrary  target  obtained  from  [1]  that  is  used  to  ensure  the  model  is 
accurate  and  the  results  are  in  agreement  with  the  reference.  This  case  is 
included  here  because  it  represents  a  complex  target  and  gives  rise  to  a  more 
complicated  returned  power  signal,  illustrated  in  the  figure  below. 


Pcol 
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Pcross 
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Pcross 


ellip  (deg)  -50  -100  tilt  (deg) 


tilt  (deg) 


Figure  IV-13  Backscattered  Power  for  an  Arbitrary  Target 


C.  FURTHER  WORK 

Even  though  most  existing  tactical  radar  systems  are  not  polarimetric, 
future  tactical  radars  are  expected  to  exploit  polarization  information.  The  current 
thesis  serves  as  a  review  of  polarimetric  investigations  as  it  provides  a 

comprehensive  summary  of  the  current  efforts  and  development  in  polarimetry 
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theories.  In  particular,  two  of  the  more  useful  and  interesting  theories  on  TD 
theorems  and  optimization  have  been  examined  at  a  level  sufficient  to  verify 
them.  Given  the  fact  that  most  of  these  polarimetry  theories  are  recent 
developments  (dating  back,  in  some  cases,  to  research  programs  begun  in  the 
late  70s),  many  aspects  of  the  proposed  theories  can  be  further  investigated, 
developed  or  verified. 

For  the  current  investigation,  a  more  complex  foliage  clutter  model  can  be 
introduced.  For  cultured  vegetation,  there  exists  some  degree  of  correlation  such 
that  the  scattering  mechanism  and  the  orientation  of  the  leaves  may  not  be 
completely  random.  The  correlation  may  be  introduced  in  the  a  angle  of  the  TD 
theorems  from  the  ‘ P'  intensities  in  the  entropy.  A  more  realistic  model  can  also 
be  achieved  if  it  is  extended  to  the  third  dimension  by  introducing  tree  trunks 
beneath  the  foliage.  The  trunks  are  expected  to  have  a  single  or  a  slightly  varying 
orientation.  Finally,  since  the  current  model  is  based  on  theoretical  analysis, 
some  realism  can  be  injected  if  it  can  be  compared  to  real  or  existing  data. 

If  existing  full  polarimetric  data  sets  can  be  obtained,  an  interesting 
problem  would  be  to  analyze  the  returns,  decompose  the  scattering  mechanisms 
and  create  a  color  plot  based  on  different  scattering  mechanisms.  An  example 
would  be  the  Pauli  decomposition  plot  given  in  Figure  111-6. 

Another  particularly  useful  implementation  of  radar  polarimetry  is  in  the 
interferometry  of  Synthetic  Aperture  Radar  (POLINSAR).  This  is  an  ongoing 
research  topic  by  current  researchers  such  as  Cloude  and  Boerner.  The  idea  is 
that  the  coherency  of  interferometry  data  sets  from  multipass  SAR  (both  temporal 
and  spatial  decorrelation)  depends  strongly  on  the  polarization  states.  Cloude 
developed  a  methodology  whereby  polarimetric  transformation  bases  can  be 
used  to  analyze  SAR  interferograms  so  that  the  coherency  can  be  optimized. 
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APPENDIX 


MATLAB  SOURCE  CODES 


SOURCE  CODE  FOR  TD  THEOREMS 


A. 

clear  all 

%Suite  of  EM-waves  (From  Ref  [12]) 
%E=[1/sqrt(2);1/sqrt(2)]; 

E=[1;0]; 

%E=[0;  1  ]; 

%E=[1/sqrt(2);(1/sqrt(2))*j] 

%E=[1/sqrt(2);-(j*1/sqrt(2))] 


%linearly  polarized  at  45degree 
%horizontally  polarized 
%verticallly  polarized 
%Left  circular  polarized 
%Right  circular  polarized 


%  Defined  alpha  angle:  scattering  mechanism 

alphaj  =  45;  %45  degrees  for  dipole 

alpha_l_rad  =  alpha_l/180*pi;  %convert  to  radians 


%Assign  random  target  location 
target_x  =round(rand*50); 
target_y  =round(rand*50); 


%Create  scene 
for  x_space=1 :50; 
for  y_space=1 :50; 


%Create  clutter  model  for  each  locations 
beta_l_rad  =  rand*pi; 

%Assign  leaves  orientation  (From  Ref  [4]) 

k_l=[cos(alpha_l_rad);sin(alphaJ_rad)*cos(betaJ_rad);sin(alpha_l_rad)*sin(beta 
_l_rad)];  %  leaves  model 


s_hh=(k_l(1,1)+  k_l(2,1))/2; 
s_vv=(k_l(1 , 1 )-  k_l(2,1))/2; 
s_x= k_l  (3,1  )/2 ; 

S(x_space,y_space,1  )=[s_hh]; 
S(x_space,y_space,2)=[s_x]; 
S(x_space,y_space,3)=[s_x]; 
S(x_space,y_space,4)=[s_vv]; 


%cl utter  scattering  matrix 
%clutter  scattering  matrix 
%cl utter  scattering  matrix 
%clutter  scattering  matrix 
%cl utter  scattering  matrix 
%cl utter  scattering  matrix 
%cl utter  scattering  matrix 


%EM  wave  depolarized  by  Clutter 
s_c  =[s_hh,s_x;s_x,s_vv]; 
e_clutter  =  s_c*E; 


%Add  effect  of  dihedral/helix/flat  plate  at  predetermined  location 
if  (x_space==target_x)  &  (y_space==target_y) 

%s_t  =[sq rt(2 )/2 , sq rt(2 )/2 ; sq rt(2 )/2 ,-sq rt(2 )/2] ;  %dihedral  at  22.5  deg 

s_t  =  [0,1  ;1 ,0];  %dihedral  at  45  deg 

%s_t  =  [0.5,-j*0.5;-0.5,-j*0.5];  %metallic  helix 
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%flat  plate 

%dihedral  at  0  degree 


%s_t=  [1 ,0;0, 1  ]; 
%s_t=  [1 ,0;0,-1  ]; 
E_t  =s_t*E; 


%Sum  clutter  and  dihedral  at  where  the  target  is  located 
E_clutter(x_space,y_space,1)=  E_t(1,1);  %+  e_clutter(1,1); 
E_clutter(x_space,y_space,2)=  E_t(2,1);  %+  e_clutter(2,1); 

else 

E_clutter(x_space,y_space,1  )=e_clutter(1 ,1 ); 
E_clutter(x_space,y_space,2)=e_clutter(2,1); 
end 


end 

end 


%Read  out  of  results 

mean(mean(S(:,:,1))) 

mean(mean(S(:,:,2))) 

mean(mean(S(:,:,3))) 

mean(mean(S(:,:,4))) 

mean(mean(E_clutter(:,:,1 ))) 

mean(mean(E_clutter(:,:,2))) 


%scattering  matrix 
%scattering  matrix 
%scattering  matrix 
%scattering  matrix 
%Back  scattered  E-field 
%Back  scattered  E-field 


%plots  for  real  values 
[x,y]  =  meshgrid([1:50]); 

Z_h  =  E_clutter(:,:,1); 

subplot(2,2,1),  surf(x,y,real(Z_h)),  title('Horizontal  Scattered  E-field') 
colorbar 


Z_v  =  E_clutter(:,:,2); 

subplot(2,2,2),  surf(x,y,real(Z_v)),  title('Vertical  Scattered  E-field') 
colorbar 


%imaginary  plots 
[x,y]  =  meshgrid([1:50]); 

Z_h  =  E_clutter(:,:,1); 

subplot(2,2,3),  surf(x,y,imag(Z_h)),  title('lmaginary  Horizontal  Scattered  E-field') 
colorbar 


Z_v  =  E_clutter(:,:,2); 

subplot(2,2,4),  surf(x,y,imag(Z_v)),  title('lmaginary  Vertical  Scattered  E-field') 
colorbar 
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B. 


SOURCE  CODE  FOR  POLARIZATION  OPTIMIZATION 


clear  all 
t=1; 

%Target  Scattering  Matrices 
S=[2j,0.5;0.5,-j]; 

%S=[1 ,0;0,1  ]; 

%S=[1 ,0;0,-1  ]; 

%S  =  [0,1;  1,0]; 

%S  =  [0.5,-j*0.5;-0.5,-j*0.5]; 

%Transitional  Matrices  (From  Ref  [1]) 
C=0.5*[1 ,0,0,0;0,1 ,0,0;0, 0,1 ,0;0,0, 0,-1]; 
X=0.5*[1 ,0,0,0;0,-1 ,0,0;0,0,-1 ,0;0, 0,0,1]; 
A=[1 ,0,0,1  ;1 ,0, 0,-1  ;0,1,1,0;0,j,-j,0]; 


%counter 


%  Arbitrary  matrix  from  Ref  [1] 

%Sphere 

%Dihedral 

%dihedral  at  45  deg 

%metallic  helix 


%Produces  the  Co-pol  Mueller 
%Produces  the  Crosspol  Mueller 
%Kronecker  expansion  matrix 


for  tilt=-0:2: 1 80; 

tilt_rad=(tilt/180)*pi; 

e=1; 

for  ellip=-45:2:45; 

ellip_rad=(ellip/180)*pi; 


%range  of  tilt  angle 
%convert  tilt  to  radians 
%counter 

%range  of  ellipticity  angle 
%convert  to  radian 


%Contruct  the  incident  E-field  based  on  the  current  tilt  and  ellipticity  (From  Ref 

[8]) 

%To  be  used  to  generate  the  Stokes  parameters 
E_i=[cos(tilt_rad)*cos(ellip_rad)- 

sin(tilt_rad)*(j*sin(ellip_rad));sin(tilt_rad)*cos(ellip_rad)+cos(tilt_rad)*(j*sin(ellip_ra 

d))]; 


%Contruct  Mueller  matrices  from  the  current  scattering  matrix 
%To  be  used  for  transforming  incident  Stokes  parameters  to  backscattered 
ones 

M  =A*(kron(S,conj(S)))*inv(A);  %general  Mueller 

Mc=C*M;  %Col-pol  matrix 

Mx=X*M;  %Cross-pol  matrix 


%Contruct  the  Stoke  Parameters  from  the  incident  E-field  (From  Ref  [11]) 

g0=abs(E_i(1  ,:)A2)+  abs(E_i(2,:)A2); 

g  1  =abs(E_i(1 , : )A2)-  abs(E_i(2, : )A2); 

g2=(abs(E_i(1  ,:)A2)-  abs(E_i(2,:)A2))*tan(2*tilt_rad); 

g3=(abs(E_i(1  ,:)A2)+  abs(E_i(2,:)A2))*sin(2*ellip_rad); 

gt=[g0;g1;g2;g3]; 

%Derived  the  backscattered  Stokes  using  the  Mueller  on  the  incident  Stokes 
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%Deriving  the  cross  pol  and  copol  powers 

power_col=transpose(gt)*Mc*gt/norm(gt,'fro'); 

power_cross=transpose(gt)*Mx*gt/norm(gt,'fro'); 

Pc(e,t)=power_col;  %data  stored  as  array 

Px(e,t)=power_cross;  %data  stored  as  array 

e=e+1;  %counter 

end 

t=t+1;  %counter 

end 

%Plots 

a=-0:2:180; 

b=-45:2:45; 

[X,Y]  =  meshgrid([-0:2: 180], [-45:2:45]); 

subplot(2,2,1),surf(X,Y,Pc),  title('Pcor),xlabel('tilt  (deg)'),ylabel('ellip  (deg)') 
subplot(2,2,2),imagesc(a,b,Pc),  title('Pcor),xlabel('tilt  (deg)'),ylabel('ellip  (deg)') 
colorbar 

subplot(2,2,3),surf(X,Y,Px),  title('Pcross'),xlabel('tilt  (deg)'),ylabel('ellip  (deg)') 
subplot(2,2,4),imagesc(a,b,Px),  title('Pcross'),xlabel('tilt  (deg)'),ylabel('ellip  (deg)') 
colorbar 

%  Readout  of  values 

[A,B]=find(Px==(max(max(Px))))  %find  ellip  &  tilt  where 

maxima  occurs 
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